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ABSTRACT 

The evolution of binaries consisting of evolved main sequence stars (1 < M^j M Q < 3.5) with white 
dwarf companions (0.7 < M w d/ M© < 1.2) is investigated through the thermal mass transfer phase. 
Taking into account the stabilizing effect of a strong, optically thick wind from the accreting white 
dwarf surface, we have explored the formation of several evolutionary groups of systems for progenitors 
with initial orbital periods of 1 and 2 days. The numerical results show that CO white dwarfs can 
accrete sufficient mass to evolve to a Type la supernova and ONeMg white dwarfs can be built up to 
undergo accretion induced collapse for donors more massive than about 2 M Q . For donors less massive 
than ~ 2 M Q the system can evolve to form a He and CO or ONeMg white dwarf pair. In addition, 
sufficient helium can be accumulated (~ 0.1 M Q ) in systems characterized by 1.6 < M<j/ M Q < 1.9 
and 0.8 < M w d/ M < 1 such that sub Chandrasekhar mass models for Type la supernovae, involving 
off center helium ignition, are possible for progenitor systems evolving via the Case A mass transfer 
phase. For systems characterized by mass ratios > 3 the system likely merges as a result of the 
occurrence of a delayed dynamical mass transfer instability. A semi-analytical model is developed 
to delineate these phases which can be easily incorporated in population synthesis studies of these 
systems. 

Subject headings: binary: close - stars: evolution - stars: mass loss - stars: cataclysmic variables - 
supernovae: general 



1. INTRODUCTION 

Close binary systems consisting of main sequence-like 
stars with white dwarf companions have long been rec- 
ognized as an important stage of evolution for under- 
standing the formation of several diverse classes of ob- 
jects. Such systems are products of an evolution of a 
main sequence-like star with a red giant or asymptotic 
giant branch companion in which significant mass and or- 
bital angular momentum have been lost. The extremely 
non conservative evolution facilitated the transformation 
of the system from a wide orbit with an orbital period of 
~ 1 year to a narrow one with a period of several days 
via a common envelope phase where the orbital energy 
released in the spiralling in process is sufficient to eject 
the common envelope (for reviews see Iben & Livio 1993; 
Taam & Sandquist 2000). For the orbital periods of post 
common envelope systems < 0.5 day, angular momentum 
losses by magnetic braking are effective (see Pylyser & 
Savonije 1988) in shrinking the orbit further, producing 
cataclysmic variable systems. On the other hand, for pe- 
riods > 0.5 day, the main sequence-like stars with masses 
> 1 M can evolve to the mass transfer stage as a result 
of envelope expansion induced by the nuclear burning in 
the main sequence star's interior. The future evolution 
of the nuclear evolved systems with donors in the mass 
range of 1 — 3 M may significantly contribute to the 
formation channels for the production of supersoft X-ray 
sources (massive white dwarfs accreting at rates suffi- 
cient for steady hydrogen burning; van den Heuvel et al. 
1992), of a class of Type la supernovae models (Li & van 
den Heuvel 1997), of ultra short period (P < 30 min) 
interacting double white dwarf AM CVn systems (Pod- 
siadlowski, Han, & Rappaport 2003), and of detached 
double white dwarf systems (Nelemans et al. 2001). 



Common to the evolution of these systems is the occur- 
rence of a phase of mass transfer on a thermal timescale. 
Until recently, this population of systems was gener- 
ally neglected in population synthesis studies since it 
was assumed that the accreting white dwarf would ex- 
pand to red giant dimensions as a result of reactiva- 
tion of hydrogen burning at high mass accretion rates 
(> 10~ 7 M yr _1 ). This expansion was hypothesized 
to lead to the formation of a second common envelope 
phase and to the eventual formation of a double de- 
generate system. However, it was pointed out by Kato 
& Hachisu (1994) that an optically thick wind can be 
driven from the surface of white dwarfs more massive 
than about 0.5 M Q , thereby stabilizing the mass trans- 
fer in the system and making new binary evolutionary 
channels possible. In this case, the photosphere of the ac- 
creting white dwarf lies within its corresponding Roche 
lobe, and the system can be stabilized preventing evo- 
lution into the common envelope phase. The existence 
of such solutions was made possible by a strong peak in 
the OPAL opacities at temperatures of about 1.6 x 10 5 K 
(Iglesias, Rogers, & Wilson 1987, 1990; Iglesias & Rogers 
1991, 1993; Rogers & Iglesias 1992) achieved in the en- 
velope of white dwarfs accreting at high rates. Recently, 
observational evidence for the accretion wind picture has 
been suggested by Hachisu & Kato (2003) based on the 
long term variability of the light curve of the transient 
supersoft X-ray source RX J0513. 9-6951 (see Reinsch et 
al. 2000). 

Our study is, in part, similar to that conducted by 
Li & van den Heuvel (1997), but differs in that we map 
out the boundaries delineating the evolutionary channels 
leading to the formation of double degenerate dwarfs, 
neutron stars by accretion induced collapse, and near 
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Chandrasekhar and sub-Chandrasekhar Type la super- 
nova models from such progenitor binary systems. As 
such, we have carried out detailed binary evolutionary 
calculations of the Roche lobe filling donor at a stage 
of evolution between the main sequence and the base of 
the giant branch at orbital periods for which nuclear evo- 
lution rather than angular momentum losses dominate. 
The evolution of the system is carried through the ther- 
mal timescale mass transfer phase to determine its ulti- 
mate outcome. To obtain a clear picture of the possible 
evolutionary histories, the mass and evolutionary state 
of the donor star and the mass of the white dwarf are 
systematically varied for systems characterized by initial 
orbital periods of 1 and 2 days. The assumptions and 
input physics underlying our calculations are described 
in §2. The detailed binary sequences and numerical re- 
sults are described and compared with a semi-analytical 
picture for the boundaries of the various evolutionary 
channels in §3. Finally, we summarize and discuss the 
implications of our results in the concluding section. 

2. FORMULATION 

The binary evolutionary sequences calculated in this 
investigation are based on a stellar evolution code de- 
veloped by Kippenhahn, Weigert, & Hofmeister (1967) 
and updated as described in Podsiadlowski, Rappaport, 
& Pfahl (2002). The stellar models are computed us- 
ing a reaction network with rates taken from Rauscher 
& Thielemann (2000, 2001) (see also Thielemann, Tru- 
ran, & Arnould 1986) and using OPAL opacities (Rogers 
& Iglesias 1992), supplemented with opacities at low 
temperatures (Alexander & Ferguson 1994), for a solar 
metallicity (Z = 0.02). 

2.1. Mass transfer 

During phases of the evolution when the donor filled its 
Roche lobe, the mass transfer rate, M tr , was calculated 
in an implicit manner. In this case, M tr is found such 
that the radius of the donor is equal to its Roche lobe 
radius, Rl approximated as (see Eggleton 1983) 
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where q is the ratio the mass of the donor, Mj, to the 
mass of the white dwarf, M w d, and A is the orbital sep- 
aration. We consider the radius-mass exponents of the 
Roche lobe £rl = d\nR^/d\nM and of the star itself 
£ = d In R/d In M in our solution method. The response 
of the Roche lobe to the mass transfer (MT) is solely 
a function of the mass ratio, whereas the response of 
the stellar radius to the MT is the function of the MT 
rate. For a given model, we tabulate values of ((M). By 
equating the Roche lobe radius and the predicted stellar 
radius, one can find the MT rate for further models. The 
MT rate solution is not necessarily unique since it can be 
multi-valued. In this case, we accept the smallest value 
of the MT rate. While the star evolves, the response of 
the stellar radius evolves as well, and we recalculate the 
table of C(M) if the predicted stellar radius differs from 
the calculated one by S\nR = 10 -4 . 

2.2. Growth of the white dwarf 



The ultimate fate of the material transferred to the 
white dwarf depends on the mass and composition of 
the white dwarf and the rate of mass transfer. In this 
study, we consider CO white dwarfs with initial masses 
between 0.7 M© and 1.15 M© and ONeMg white dwarfs 
more massive than 1.15 M Q . The actual growth of the 
white dwarf is determined by the amount of accreted 
hydrogen which is eventually converted to helium and 
heavier elements. In the case of rapid mass transfer, the 
white dwarf accumulates matter at the rate 

M„ ~ 7.5 x 10~ 7 (Af wd / M© - 0.4) M© yr" 1 , (2) 
with the remaining matter radiatively driven away in a 
strong optically thick superwind from the white dwarf 
surface (Hachisu, Kato, & Nomoto 1999). The accu- 
mulation ratio, tjh = M cr /M tr < 1 . For transfer 
rates less than M cr , but greater than about M\ ow ~ 
10 -7 M© yr _1 , hydrogen burns steadily to helium {t]h ~ 
1). If the mass transfer rate varies in the range between 
~ 3 x 10~ 8 M© yr _1 and M\ ow , the white dwarf experi- 
ences mild recurrent hydrogen shell flashes. For the mass 
transfer rates below ~3x 10~ 8 M© yr _1 , strong unsta- 
ble hydrogen shell flashes result in nova explosions, and 
we assume that the evolution is fully non conservative 
with all the transferred material ejected from the sys- 
tem. In this regime, the white dwarf mass is constant 
(i.e., t]h = 0). 

In our calculations, r\n is based on the work of Pri- 
alnik & Kovetz (1995) who considered the accretion of 
hydrogen-rich matter onto white dwarfs of varying mass 
and central temperature. Since the accumulation ratios 
do not sensitively depend on the thermal state of the 
white dwarf, fits of their results were used for accre- 
tion rates ranging from 10 _9 M©yr _1 to 10 _6 M©yr _1 for 
their hot white dwarf models. 

In the mass accretion rate regime where hydrogen-rich 
matter is retained, helium will naturally be accumulated. 
For a sufficiently massive helium layer the helium will be 
ignited. The critical mass necessary for helium ignition is 
dependent on the rate at which helium is processed in the 
hydrogen burning shell Mh c (Kato & Hachisu 1999). For 
rates less than 1.3 x 10~ 6 M© yr~ x , helium shell burning 
is unstable in the white dwarf envelope. The accumula- 
tion ratio, defined as the ratio of the processed material 
remaining after one cycle of a helium flash to the ignition 
mass, has been estimated by Kato & Hachisu (1999) as 

-5.9<logM Hc <-5 , * 
^ -0.175(logM H c + 5.35) 2 + 1.05), (6) 
-7.3 < logM He < -5.9 

where Mhc is units of Af© yr^ 1 . We note that the He 
shell is always unstable whenever a significant amount of 
hydrogen is processed to helium. 

If Mae falls below the lower range in Eq. |31 the helium 
layer is not sufficiently hot to ignite helium shell burning 
in a thin mass layer. In this case, helium is accumulated 
in the layer until the increasing density leads to the he- 
lium ignition in a thick mass layer of about 0.1 M©. 
This ignition of the helium layer likely leads to the igni- 
tion of the CO or ONeMg core disrupting the star as a 
sub-Chandrasekhar mass Type la supernova (see Taam 
1980; Livne & Glasner 1991; Woosley & Weaver 1994; 
Garcia-Senz, Bravo, & Woosley 1999). 
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2.3. Orbital evolution 

The evolution of the orbital separation is dependent on 
the prescription for mass and angular momentum loss. 
Since we consider systems characterized by orbital peri- 
ods (P > 1 day) exceeding the bifurcation period (Pyl- 
yser & Savonije 1988) the only loss of orbital angular 
momentum from the system reflects that carried by the 
ejected matter. In this study we assume that the mass 
lost in the radiatively driven wind carries the specific or- 
bital angular momentum of the white dwarf. Hence, the 
rate of angular momentum loss is given by 



M d 
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, r -M wind (4) 
M wd M d + M wd 

where the orbital angular momentum, J, of the system 
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and M W ind is the mass loss rate from the system. 

The orbital and stellar evolution coupled to the stellar 
response to mass loss implicitly provides the constraint 
imposed on the system for determining the rate of mass 
transfer as described above. 

3. RESULTS 

In order to facilitate an understanding of the detailed 
numerical results we, first, approximately estimate the 
boundaries delineating the possible outcomes of the evo- 
lution in the next subsection. A description of the de- 
tailed binary evolutionary sequences follows in the sub- 
sequent subsections. 

3.1. Semi- analytical picture 

The evolutionary fate of a given initial binary system 
is critically dependent on the rate of mass transfer during 
the Roche lobe overflow phase. As a result, an estimate 
for the average mass transfer rate is required to deter- 
mine the approximate boundaries separating the possi- 
ble outcomes. In the following, we take the average mass 

transfer rate as 

M = AM/t th , (6) 

where tth is the thermal time-scale, on which the mass 
transfer operates, and AM is the mass lost by the donor 
during this mass transfer event. The thermal mass trans- 
fer from a MS star or a giant star takes place when the ac- 
cretor is more massive star and terminates when M d (f) = 
M d -AM = M wd +J7 H (M)?7He(MHe)AM. In this descrip- 
tion, only AM = (M d - M wd )/(1 + 7? H (M)r7Hc(M H o)?7Ho) 
can be lost during the thermal mass transfer phase. In 
the case of the star evolving through the Hertzsprung 
Gap (HG), the thermal time mass transfer is driven also 
by the star's own expansion, and complete envelope can 
be lost (AM = M env ) during this phase, however we do 
not treat it separately in our simplified picture. 

The thermal timescale tth on which AM is lost, is ap- 
proximately given as 



, _ , M d /Mq AM/ Mq 
tth - tth,© Rd/Re Ld/L(D 

_ , (M d /M )- 2 - 25 AM/M 
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HG, since they maintain the same luminosity as at the 
end of the MS), tth,© is the thermal time-scale of the Sun. 
In this simplified estimate, we do not follow the evolution 
of the star through the MT, and therefore use only the 
value of tth at the onset of the MT. The stellar radius, 
R, is identified then with the radius of the Roche lobe 
and is a function of the orbital period of the binary and 
the mass ratio of the system. This yields a very simple 
description for the average mass transfer rate as 



M = k 



R d / R (M d /M ) 

tth,© 



2.25 



(8) 



Here we used the mass-luminosity relation for main se- 
quence stars, L oc M 3,25 , (this is also correct for stars at 



where k is some factor of the order of unity, which can 
be taken from the detailed calculations. 
The final mass of the accretor is M wd j = M wd + 

r)H_(M)r)n c (Mii c )Mtth (for the complete description for 
?7h and rjn e see §2.2). This simple approximate approach 
permits an exploration of the parameter space for the 
final product of the binary evolution as a function of ini- 
tial donor mass, white dwarf mass, and orbital period. 
As representative examples we illustrate the boundaries 
for the possible evolutionary fates of systems character- 
ized by initial orbital periods of 1 and 2 days in Figs. 1 
and 2 respectively. For P = 1 day, where the thermal 
mass transfer occurs near the main sequence we adopt 
tth,© ~ 3 x 10 7 yr. On the other hand, a 1 M Q star in 
the HG has a thermal timescale of 1.8 x 10 7 yr. This 
corresponds to tth,© ~ 6 x 10 7 yr in equation (7). For 
both Fig. 1 and 2 a value of k = 1.5 was used based on 
the results of our detailed calculations 

The parameter range for which a white dwarf can ac- 
cumulate sufficient mass to either evolve to a 1.4 M Q 
white dwarf required for a Type la supernova event (SN 
la) or an accretion induced collapse (AIC) is enclosed 
within the solid boundaries. The initial mass of the 
white dwarf is the distinguishing characteristic in these 
evolutionary channels since CO white dwarfs with ini- 
tial masses < 1.15 M Q explode upon the ignition of car- 
bon, whereas ONeMg white dwarfs with initial masses 
> 1.15 M Q collapse to form neutron stars as a result of 
electron capture processes. The mass and composition 
of the initial white dwarf component of the system is 
necessary, but not sufficient to ensure these evolution- 
ary fates, since the donor must supply matter at a rate 
such that the white dwarf can grow significantly. Only 
donors of an intermediate mass range can satisfy such 
a condition since the evolution is truncated for massive 
donors by the onset of a delayed dynamical mass trans- 
fer instability (Webbink 1985; Hjellming 1989) and for 
low masses by the occurrence of the nova phenomenon. 
Specifically, the upper limit to the donor's mass is a con- 
sequence of the change in the donor's response when the 
steep entropy profile of the outer envelope layers is re- 
moved. In this case stellar contraction is replaced by 
stellar expansion as the flat entropy profile of the inte- 
rior layers is exposed when sufficient matter is removed 
on a thermal timescale. Since the Roche lobe is contract- 
ing during this phase, the mass transfer process becomes 
dynamically unstable. Based on the results of our de- 
tailed calculations wc adopt a critical mass ratio of 2.9 
(for P = 1 day) and 3.1 (for P = 2 days) above which the 
delayed dynamical instability (DDI) is found. For sys- 
tems which enter the DDI phase, it is highly likely that 
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Fig. 1. — Possible evolutionary fates of systems characterized by 
donors of initial mass M d and white dwarfs of initial mass M wd 
in units of Mq for an orbital period of 1 day based on a semi- 
analytical model (sec text). The solid curve delineates the area 
where the white dwarf can accrete sufficient mass for near Chan- 
drasekhar mass models of SN la and AIC models of neutron star 
formation. The dashed curve represents the condition of the oc- 
currence of a DDI (mass ratio q=3.1). Double white dwarf systems 
are formed at locations marked DWD. For reference, donors which 
undergo mass transfer above the upper dotted line lie close to the 
main sequence, those below are in the Hertzsprung gap, and those 
below the lower dotted line are red giants. Above the short dashed 
- long dashed line the mass transfer leads to the development of a 
superwind from the white dwarf. For stars that initiate mass trans- 
fer near the main sequence, but for which the rate is insufficient to 
lead to the superwind or to the SN la phase, a sub-Chandrasekhar 
SN is a possible outcome. 
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Fig. 2. — Possible evolutionary fates of systems characterized by 
donors of initial mass M d and white dwarfs of initial mass M wd 
in units of Mq for an orbital period of 2 days based on a semi- 
analytical model (see text). In addition to the boundaries described 
in Fig. 1, a lower dashed curve separates systems which enter 
into the common envelope phase (corresponding to mass ratios 
> 1.3) from those systems which evolve to form double white dwarf 
systems without requiring a common envelope phase. 



the system will merge, leaving behind a rapidly rotating 
remnant. This upper limit, based on the DDI criterion, 
overestimates the donor's mass for systems characterized 
by white dwarfs less massive than 1.05 M Q since donors 
in the mass range between 1.8 M Q and 3.3 M either 
are not sufficiently massive to build up the white dwarf 
to 1.4 M (even assuming an accumulation efficiency of 
unity) or because the efficiency of accretion is strongly 
reduced by the radiatively driven wind. A lower limit 
on the donor mass for the SN la and AIC evolutionary 
scenario is determined by the occurrence of strong un- 
stable hydrogen shell flashes in the envelope of the white 
dwarf. In particular, the average mass transfer rate as 
determined by the thermal timcscalc of the donor is insuf- 
ficiently high for stars < 1.7 M Q to prevent the ejection 
of significant mass from the system via nova explosions. 

The remaining combinations of donor and white dwarf 
masses primarily lead to the formation of a double white 
dwarf (DWD) systems. In these cases the mass transfer 
rate is insufficient to lead to the growth to the SN la or 
AIC phase. We note that there exists a narrow range 
of parameter space (Md ~ 1.7 Mq and M wd - ~ 0.8 — 
1 Mq) for which a significant layer of helium may be 
accumulated, giving rise to an evolutionary channel in 
which a sub-Chandrasekhar mass supernova model may 
be viable. 

The range in donor and white dwarf masses delineating 
the evolutionary fates of systems in Fig. 2 for an initial 
orbital period of 2 days at the onset of mass transfer are 
similar to those described in Fig. 1, allowing for the pos- 
sibility that white dwarfs can be significantly built up. 
For this greater orbital period, the main sequence like 
donors are evolved to a greater extent and the onset of 
mass transfer occurs after the formation of a helium core 
when the donor is in the Hertzsprung gap or at the base 
of the giant branch. As is evident from Fig. 2 there ex- 
ists a region where the system can enter into a common 
envelope phase as a result of evolution to the giant stage. 
The lower dashed line denotes this region corresponding 
to the condition that the mass ratio equals 1.3 as based 
on our detailed numerical calculations. For mass ratios 
greater than 1.3 the system enters into the common enve- 
lope phase. In this case the donor has a well defined core- 
envelope structure typical of red giant stars. A double 
degenerate system consisting of a He white dwarf with a 
CO or ONeMg white dwarf in a short period orbit (P < 1 
hour) may result depending on the particular evolution- 
ary state of the donor at onset of the common envelope 
stage (Sandquist, Taam, & Burkert 2000). An additional 
difference between the systems depicted in Fig. 1 and in 
Fig. 2 is the absence of systems which can evolve to a 
sub-Chandrasekhar mass SN. This reflects the fact that 
the regime where a sufficient helium mass layer accumu- 
lates on the white dwarf occurs for systems in which the 
mass transfer takes place only when the donor is close 
to the main sequence. For donors of more advanced evo- 
lutionary stages in mass transferring systems at initial 
orbital periods of 2 days, the timcscalc for evolution is 
so short that significant accumulation of helium required 
for a sub Chandrasekhar model is not found. 

3.2. Detailed calculations 

We have computed the evolution of binary systems ini- 
tially consisting of main sequence-like stars of masses 1 - 
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Fig. 3. — The evolutionary fates of the binary sequences cal- 
culated for an initial orbital period of 1 day and donors of initial 
mass M d and white dwarfs of initial mass M wd in units of Mq. 
The symbols representing the fate of the system as a double white 
dwarf (DWD), a near Chandrasckhar mass supernova model (SN 
la), a helium ignition sub Chandrasckhar model (He ignition), ac- 
cretion induced collapse (AIC), and an evolution into the common 
envelope stage via a delayed dynamical mass transfer instability 
(CE DDI) are shown. 
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Fig. 4. — The evolutionary fates of the binary sequences calcu- 
lated for an initial orbital period of 2 days and donors of initial 
mass M d and white dwarfs of initial mass M wd in units of Mq. 
The symbols representing the fate of the system as a double white 
dwarf (DWD), a near Chandrasckhar mass supernova model (SN 
la), accretion induced collapse (AIC), and an evolution into the 
common envelope stage via a dynamical instability (CE) leading 
to merger and to possibly a DWD system (CE — * DWD), and 
delayed dynamical mass transfer instability (CE DDI) arc shown. 



3.8 M© and white dwarfs of masses 0.7—1.2 M© through 
the thermal mass transfer phase. The evolutionary state 
of the donor stars was chosen such that the Roche lobe 
overflow phase was initiated at either an orbital period of 
1 or 2 days. To sample the parameter space adequately, 
a total of 65 evolutionary sequences were calculated with 
28 and 37 sequences calculated for systems at initial or- 
bital periods of 1 and 2 days respectively. The systems 
were evolved until the accreting white dwarf had reached 
1.4 M©, the donor had evolved to a point where it was 
clear that the system would evolve into a DWD system, 
or the system enters into a common envelope phase as a 
result of either a dynamical or delayed dynamical mass 
transfer instability. 

A visual summary of the fate of all the calculated evo- 
lutionary sequences is presented in Figs. 3 and 4 for 
initial orbital periods of 1 and 2 days respectively. Upon 
inspection of the boundaries displayed in Figs. 1 - 4, it 
is clear that the results based on the semi-analytical ap- 
proach compare very favorably with the detailed binary 
evolutionary computations. The quantitative results of 
representative evolutionary sequences are listed in Table 
1 and 2 for initial orbital periods of 1 and 2 days respec- 
tively. Here, the initial mass of the white dwarf, M w d ; o, 
mass of the white dwarf at the end of the calculations, 
M w d-j, the initial mass of the donor, M^.q, mass of the 
donor at the end of the calculations, Md-j, the minimum 
orbital period that the binary system has reached during 
the computed evolution, P m i n , the orbital period at the 
end of the calculations, Pf, the evolution time from the 
onset of the mass transfer phase, t tr , and the average 
mass transfer rate, M tr are listed. 

In the following, we divide the description of the de- 
tailed numerical results into those sequences for which 



growth of the accreting white dwarf is sufficient to lead to 
a near Chandrasekhar or sub-Chandrasekhar mass model 
for SN la and to an AIC phase and sequences leading 
to the formation of double degenerate systems or to a 
merged object. 

3.2.1. Growth of white dwarf to SN la or AIC phase 

The progenitor systems which favor the growth of the 
white dwarf to a near Chandrasekhar mass are charac- 
terized by massive white dwarfs (M w d ; o ?J 0.8 M©). It 
can be inferred from Tables 1 and 2 that the efficiency of 
mass accretion by the white dwarf, as defined by the ratio 
of the mass accreted by the white dwarf to the mass lost 
by the donor, is a strong function of the mass and evolu- 
tionary state of the donor. Specifically, the efficiency of 
white dwarf growth ranges from about 10% to 70% over 
the calculated grid with the highest efficiencies obtained 
for donors with masses ~ 2 M© . Clearly the mass trans- 
fer rates for these donors are in the regime where steady 
hydrogen burning takes place, but yet the superwind is 
not so effective in driving a significant amount of mat- 
ter from the system. That is, the thermal timescale mass 
transfer associated with donors of ~ 2 M Q favors a higher 
efficiency of material accumulation on the white dwarf in 
comparison to other donors. In contrast, the progenitor 
systems which lose a significant fraction of their mass via 
the radiatively driven wind from the white dwarf surface 
during the evolution generally are characterized by sys- 
tems with massive donors > 2.5 M©. In fact, the degree 
to which systems lose mass is exemplified by the systems 
which evolve to the SNIa or AIC phase. For these sys- 
tems, the total mass loss from the system can lead to the 
presence of ~ 0.4 — 1.9 M© surrounding the system at 
the time the white dwarf has been built up to 1.4 M©. 
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Fig. 5. — Evolution of a binary characterized by Md;0 = 3-2 Mq, 
M„,j.(| = 1.2 Mq at an initial orbital period of 2 days. Upper 
panel: The temporal evolution of the orbital period. Middle panel: 
The temporal evolution of the mass transfer rate. Lower panel: 
The variation of the mass of the donor and white dwarf during the 
evolution. 



An example of such an evolution characterized by an 
ONeMg white dwarf of M w d-o = 1.2 M Q and a donor 
of Md-o = 3.2 M© at an initial orbital period of 2 days 
is illustrated in Fig. 5. In this case the mass transfer 
rates rise rapidly to ~ 10~ 5 M© yr -1 within 10 5 yrs and 
averages about 6.8 x 10 -6 M© yr -1 over the entire evolu- 
tion. The efficiency of the mass accretion onto the white 
dwarf is low in this sequence amounting to ~ 10%. The 
orbital period initially decreases by more than a factor of 
2 reaching a minimum period of 0.8 days after 2 x 10 5 yrs 
before increasing to a period of 1.1 day at which point the 
white dwarf has increased to 1.4 M© and the donor has 
decreased to 1.1 M Q . The mass loss in such a sequence 
has been extensive with about 1.9 M© enveloping the 
system. It is expected that the system would undergo an 
accretion induced collapse further widening the orbit as 
the gravitational mass of the compact object is reduced 
by about 0.1 — 0.2 M Q . With the additional expansion of 
the donor, the system will evolve to longer orbital periods 
and enter the low mass X-ray binary phase. 

We have also found that systems with donors and 
white dwarfs in a narrow range (1.6 < Md-o/ M a < 1.9, 
0.8 < M w d-o/ M < 1) can lead to the accumulation of 
a sufficient layer of helium mass on the CO white dwarf 
required for the initiation of an off center helium detona- 
tion. The conditions required for the growth of the he- 
lium layer are a function of the evolutionary state of the 
donor since the strong helium flash regime was found to 
only occur for donors which undergo the MT phase near 
the main sequence at short orbital periods (P ~ 1 day). 

3.2.2. Formation of DWD or single merged object 

Although the evolutionary sequences were carried out 
only through the thermal timescale mass transfer phase, 
the outcome of an evolution as a DWD system can be 
inferred once the mass transfer rates in the system have 
decreased to the point when strong hydrogen shell flashes 
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Fig. 6. — Evolution of a binary characterized by M^o = 3.4 Mq, 
M w d o = 1 Mq at an initial orbital period of 2 days for a system 
that evolves to a delayed dynamical mass transfer instability. Up- 
per panel: The temporal evolution of the orbital period. Middle 
panel: The temporal evolution of the mass transfer rate. Lower 
panel: The variation of the mass of the donor during the evolu- 
tion. The white dwarf mass is not shown since it increases by only 
0.02 M . 



are expected to occur in the white dwarf envelope. Fig- 
ures 3 and 4 illustrate that DWD systems can form from 
progenitor systems with either a low mass (< 1.4 M©) or 
high mass (> 2.6 M©) donor. For the lower mass donors, 
the temporal evolution of the mass transfer rates lead 
to nova explosions over part of the evolution such that 
the white dwarf does not build up to 1.4 M©. On the 
other hand, there is a narrow range for high mass donors 
where, although the white dwarf can be built up signif- 
icantly, the donor becomes less massive than the white 
dwarf during the evolution such that the mass transfer 
rates have decreased to the extent that accretion is pre- 
vented by the occurrence of nova explosions. The DWD 
systems formed via this channel are expected to have 
long orbital periods [P > 1 days). 

The evolutionary channel involving a CE phase may 
also produce DWD systems. The calculations reveal that 
this evolutionary channel is restricted to the mass trans- 
fer taking place after the donor has evolved through the 
Hertzsprung gap. The response of the star with a deep 
convective envelope leads to a mass transfer instability 
and to the evolution into the CE phase. The outcome 
of this evolution is not well understood, but if the sys- 
tem survives it is likely to have a short orbital period 
(P < 1 hr). The common envelope calculations carried 
out by Sandquist, Taam, & Burkert (2000) indicate that 
evolution into the CE near the base of the giant branch 
may lead to the successful ejection of the common enve- 
lope if the degenerate helium core is more massive than 
- 0.2 - 0.25 M . 

For more massive donors, the thermal timescale mass 
transfer phase enters into a period where the mass trans- 
fer rate accelerates rapidly, indicating the onset of the de- 
layed dynamical mass transfer instability. For example, 
a sequence characterized by Md = 3.4 M©, M w d = 1 M© 
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at an initial orbital period of 2 days is illustrated in Fig. 
6. The mass transfer rate accelerates by three orders 
of magnitude after ~ 4.9 x 10 4 years with the rapid in- 
crease occurring over a timescale of < 200 years. In 
these sequences, the mass transfer initially takes place 
on a thermal timescale of the donor until the outer ra- 
diative envelope layer has been removed. Further mass 
loss leads to the exposure of the layers characterized by 
a flat entropy profile and expansion of the donor results. 
Since the Roche lobe during this phase is contracting, 
the system evolves into a CE phase. In this case, merger 
of the system is likely. The occurrence of this insta- 
bility occurs at a mass ratio of about 3.1 for donors in 
the Hertzsprung gap and about 2.9 for donors near the 
main sequence, confirming earlier estimates by Hjellming 
(1989). 

4. SUMMARY 

The fate of binary systems which undergo a phase of 
mass transfer on a thermal timescale has been investi- 
gated for binaries composed of main sequence-like donors 
with white dwarf companions. Allowing for the possibil- 
ity of an optically thick wind driven from the white dwarf 
surface gives rise to evolutionary channels in which the 
white dwarf accretes sufficient mass as required for cither 
a near Chandrasekhar mass or sub-Chandrasekhar mass 
model for Type la supernovae or the formation of neu- 
tron stars via an accretion induced collapse model. The 
range for donor and white dwarf masses delineating the 
formation of these objects from those producing a pair 
of degenerate dwarfs has been determined by detailed 
computations, and reproduced by a simple semi-analytic 
model, for progenitor systems characterized by initial or- 
bital periods of 1 and 2 days. 

Specifically, we have found that CO white dwarfs can 
accumulate sufficient matter to produce conditions ripe 
for initiation of a central carbon deflagration/detonation 
supernova explosion provided that their initial masses 
are greater than about 0.8 Mq and their donor is more 
massive than about 2 Mq. An upper limit to the donor 
mass is set by the onset of a delayed dynamical instabil- 
ity which occurs at a mass ratio of ~ 3. In addition, the 
ONeMg white dwarfs (with masses > 1.15 M©) may ac- 
crete sufficient matter to collapse to form a neutron star 
by electron capture processes for binary systems with 
donors in the same mass range. During these phases, the 
mass transfer rates are sufficiently high that hydrogen 
burning provides the bulk of the energy generation such 
that the sources are likely to be observed as supersoft X- 
ray sources (van den Heuvel et al. 1992; Kahabka & van 
den Heuvel 1997). These evolutionary results are similar 
to those reported by Li & van den Heuvel (1997) and 
extend the lower bound of the progenitor white dwarf 
masses from 0.9 M Q to about 0.8 M Q and less, confirm- 
ing similar results obtained by Langer et al. (2000) 1 . 
Although Li & van den Heuvel (1997) did not consider 
accretion induced collapse in their study, their results are 
also applicable for this neutron star formation path as 
well. In these cases, the amount of matter surrounding 

1 In a very recent study on SN la progenitors an even lower 
bound for WD masses was reported - 0.67 Mq (Han & Podsiad- 
lowski, astro-ph/0309618 ). This value, however, was based upon 
the full accretion of hydrogen rich material even for M < M\ avl . 



the system as a result of the wind loss from the white 
dwarf surface may be as high as 1 — 2 M Q . The ex- 
istence of such circumbinary material is not dissimilar 
to that expected for complementary evolutionary chan- 
nels leading to a SN la involving mass transfer from an 
asymptotic giant branch star to its massive white dwarf 
companion. We note that additional hydrogen rich mat- 
ter may also be lost from the system resulting from the 
interaction of the supernova shell with the donor (Ma- 
rietta, Burrows, & Fryxell 2000). In the past, one of 
the arguments put forth against such accreting models 
is the lack of hydrogen emission in the spectra of SN 
la. However, the recent discovery of hydrogen emission 
in the observations of the type la supernova SN2002ic 
by Hamuy et al. (2003) provides some support for such 
models. Although Hamuy et al. (2003) attribute the 
hydrogen emission from circumstellar matter surround- 
ing an asymptotic giant branch star, the wind present 
in the accreting white dwarf model in the short period 
systems investigated here suggests that sufficient matter 
surrounding the system may be a characteristic of these 
models as well. 

Our detailed systematic investigation has also led to an 
identification of a channel where a sufficient helium-rich 
mass layer accumulates (~ 0.1 Mq) below the hydrogen- 
rich layer such that an off center helium detonation may 
be ignited in the white dwarf envelope. The propagation 
of the nuclear burning front into the core region, en- 
hanced by geometrical focusing, may lead to the inciner- 
ation of the entire white dwarf (Woosley & Weaver 1994) 
for a sub-Chandrasekhar mass star. The range of masses 
of the progenitor systems as well as the evolutionary 
state of the donor which follow this evolutionary channel 
are narrow in the parameter range studied in this paper. 
That is, only donors in the mass range between 1.6 and 
1.9 Mq which transfer matter to their white dwarf com- 
panion (with masses in the range 0.8 < M w d/ M Q < 1) 
while close to the main sequence are viable. Such limited 
properties of the progenitor systems may be consistent 
with the fact that such models are expected to be mi- 
nor contributors to the SN la rate since these models are 
subluminous compared to models based on the carbon 
deflagration/detonation near Chandrasekhar mass mod- 
els. 

For donors less massive than ~ 2 M Q less matter is 
transferred to the white dwarf, leading to systems com- 
posed of a He white dwarf with a CO or ONeMg white 
dwarf companion. Of these systems, those that undergo 
a dynamical mass transfer instability evolve into the com- 
mon envelope phase. For those systems which survive, 
the systems are likely to emerge at short orbital periods 
(P < 1 h) , providing an alternative evolutionary channel 
for the formation of AM CVn binary systems in addition 
to the channel involving the stable mass transfer evolu- 
tion of evolved secondaries in cataclysmic variable sys- 
tems (see Podsiadlowski, Han, & Rappaport 2003). On 
the other hand, the systems which do not enter into the 
common envelope phase produce double degenerate sys- 
tems characterized by orbital periods which are greater 
than about 1 day. 

The systems which undergo an accretion induced col- 
lapse likely evolve to a low mass X-ray binary phase in 
which a main sequence-like star with an evolved core 
transfers mass to its newly formed neutron star com- 
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panion. Investigations of this phase have recently been 
carried out by Sutantyo & Li (2000), who considered the 
accretion induced collapse scenario, as well as the in- 
termediate mass X-ray binary scenario (also studied by 
Podsiadlowski et al. 2002). The system evolves to be- 
come a binary millisecond pulsar with a helium white 
dwarf companion in a short period system. Our results 
suggest that such systems which form via AIC process 
are likely to have orbital periods greater than about 1 
day. 

Finally, we point out that the semi-analytical model for 



Alexander, D. R., & Ferguson, J. W. 1994, ApJ, 437, 879 
Eggleton, P. P. 1983, ApJ, 268, 368 

Garcia-Senz, D., Bravo, E., & Woosley, S. E. 1999, A&A, 349, 177 

Hachisu, I., Kato, M., & Nomoto, K. 1999, ApJ, 522, 487 

Hachisu, I., & Kato, M. 2003, ApJ, 590, 445 

Hamuy, M. et al. 2003, Nature, in press 

Hjellming, M.S. 1989, Ph.D. thesis, University of Illinois 

Iben, I. Jr., & Livio, M. 1993, PASP, 105, 1373 

Iglesias, C. A., & Rogers, F. 1991, ApJ, 371, L73 

Iglesias, C. A., & Rogers, F. 1993, ApJ, 412, 752 

Iglesias, C. A., Rogers, F. J., & Wilson, B. G. 1987, ApJ, 322, L45 

Iglesias, C. A., Rogers, F. J., & Wilson, B. G. 1990, ApJ, 360, 221 

Kahabka, P., & van den Heuvel, E. P. J. 1997, ARA&A, 35, 69 

Kato, M., & Hachisu, I. 1994, ApJ, 437, 832 

Kato, M., & Hachisu, I. 1999, ApJ, 513, L41 

Kippcnhahn, R., Weigert, A., & Hofmeister, E. 1967, in Methods 

in Computational Physics, Vol. 7, ed. B. Alder, S. Fcrnbach, & 

M. Rothenberg (New York: Academic), 129 
Langer, N., Deutschmann, A., Wellstein, S., & Hoflich, P. 2000, 

A&A, 362, 1046 
Li, X. D., & van den Heuvel, E. P. J. 1997, A&A, 322, L9 
Livnc, E., & Glasner, A. S. 1991, ApJ, 370, 272 
Marietta, E., Burrows, A., & Fryxell, B. 2000, ApJS, 128, 615 
Nelemans, G. Yungelson, L. R., Portegies Zwart, S. F., & Verbunt, 

F. 2991, A&A, 365, 491 
Podsiadlowski, Ph., Han, Z., & Rappaport, S. 2003, MNRAS, 340, 

1214 



the white dwarf binaries that we have developed in this 
paper can be easily incorporated into population synthe- 
sis investigations. We plan to carry out such calculations 
in the future in order to assess the importance of the for- 
mation channels associated with the thermal timescale 
mass transfer phase in these systems. 
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Table 1. The binary parameters and outcomes of representative 
model sequences for systems in which the onset of mass transfer oc- 
curs at a orbital period of 1 day. 



M wd;0 


M d . 


M wd;f 


M d . f 


Pnin {d\ 


Pf M 


A* tr [yr] 


M tl [McJyr] 


outcome 


0.7 


2.0 


1.34 


0.89 


0.42 


0.5 


5.5-10 6 


2 • 10~ 7 


DWD 


0.8 


1.5 


0.81 


0.75 


0.82 


1.02 


7.9-10 7 


9.5 ■ 10~ 9 


DWD 


0.8 


1.7 


0.9 


1.25 


0.69 


0.69 


9.5-10 6 


4.7- 10~ 8 


subCh 


0.8 


2.0 


1.4 


1.17 


0.52 


0.55 


3.9-10 6 


2.2 ■ 10~ 7 


SNIa 


1.0 


1.0 


1.02 


0.90 


1.0 


1.02 


1.1-10 7 


9 ■ 10~ 9 


DWD 


1.0 


2.5 


1.4 


1.37 


0.51 


0.51 


1.7- 10 s 


6.6 • 10~ 7 


SNIa 


1.2 


2.0 


1.4 


1.70 


0.86 


0.86 


3.3-10 6 


9.1 • 10~ 8 


AIC 


1.2 


3.0 


1.4 


1.64 


0.5 


0.5 


2.9-10 6 


4.7- 10~ 7 


AIC 



Note. — The columns denote: the mass of the white dwarf, M wc j ; o an d of the donor M^q at the onset of the MT (in Mq); the mass of 
the white dwarf M wc j;f (in Mq), the mass of the donor Ma ; f (in Mq), and the period Pf (in days) at the end of the computation; P m in (in 
days) is the minimum period of the binary system during the computed interval of evolution, Aitr [y r ] is the time of the MT phase, and Mtr (in 
Mq yr — 1 ) is the corresponding average MT rate. The possible outcomes are: DWD: double white dwarf system, SNIa - supernova type la, subCh 
- sub-Chandraskehar supernova, AIC - accretion-induced collapse. 



Table 2. The binary parameters and outcomes of representative model 
sequences for systems in which the onset of mass transfer occurs at 
an orbital period of 2 days. 





M d . 


M wd;f 


M d:i 


Pmin [d] 


P f [d] 


Afe [yr] 


Mtr [M Q /yr] 


outcome 


0.8 


2.0 


1.4 


1.25 


1.05 


1.08 


2.-10 6 


3.8 -10" 7 


SNIa 


0.8 


2.4 


1.32 


0.32 


0.58 


4.2 


9.1-10 6 


2.3 • 10~ 7 


DWD 


1.0 


1.4 


1.0 


0.2 


2.0 


0.02 a 






CE^DWD 


1.0 


2.0 


1.4 


1.46 


1.44 


1.44 


2.6-10 6 


2.1 • 10~ 7 


SNIa 


1.0 


2.4 


1.4 


1.4 


1.1 


1.1 


8.4-10 5 


1.2 • 10~ 6 


SNIa 


1.0 


2.8 


1.4 


0.6 


0.7 


1.7 


7.7-10 5 


2.9 ■ 10~ 6 


SNIa 


1.2 


2.0 


1.4 


1.67 


1.7 


1.7 


2.5-10 6 


1.3 • 10~ 7 


AIC 


1.2 


2.6 


1.4 


1.97 


1.4 


1.4 


3.6- 10 5 


1.7 • 10~ 6 


AIC 


1.2 


3.2 


1.4 


1.1 


0.8 


1.1 


3.110 5 


6.8 • 10~ 6 


AIC 



Note. — The columns arc the same as for Tabic 1; the additional possible outcome, CE— »DWD, corresponds to systems that possibly survive 
the common envelope phase and form a short period double white dwarf system. 
a dctcrmincd by equating the orbital energy to the binding energy of the donor's envelope 



